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J5 . ABSTRACT 

I re-examine the brightness temperature problem in PKS 0405-385 which is an extreme 
04 ' intra-day variable radio quasar with an inferred brightness temperature of ~ 5 x 10^'* K 

at 5 GHz, well above the Compton catastrophe limit of ~ 10^^ K reached when the 
synchrotron photon energy density exceeds the energy density of the magnetic field. If 
one takes into account the uncertainty in the distance to the ionized clouds responsible 
for interstellar scintillation causing rapid intra-day variability in PKS 0405-385 it is 
' possible that the brightness temperature could be as low as ^ 10^^ K at 5 GHz, or even 

, lower. The radio spectrum can be fitted by optically thin emission from mono-energetic 

electrons, or an electron spectrum with a low-energy cut-off such that the critical 
I frequency of the lowest energy electrons is above the radio frequencies of interest. If 

CO . one observes optically thin emission along a long narrow emission region, the average 

' energy density in the emission region can be many orders of magnitude lower than 

fH I calculated from the observed intensity if one assumed a spherical emission region. 

Q^. I discuss the physical conditions in the emission region and find that the Compton 

I ' catastrophe can then be avoided using a reasonable Doppler factor. I also show that 

O ^ MeV to 100 GeV gamma-ray emission at observable flux levels should be expected 

-4— » ' from extreme intra-day variable sources such as PKS 0405-385. 

' Key words: galaxies: active - galaxies: jets - gamma-rays: theory - radiation mech- 

anisms: non-thermal - quasars: PKS 0405-385. 



1 INTRODUCTION convert the observed 4.8 GHz flux to intensity, and obtain a 

/TT-^trN variability brightness temperature of Tvar ~ 10^^ K. In 
Rapid variability m mtra-day variable (IDV) sources is a . , i • i . i . . 

^ ' this source, however, Kedziora-Chudczer et al. (1997) mter- 
lone-standme problem as it implies apparent brishtness tern- , . , i ^ ■ ^- i ■ ^ ^ n 

° • u J- • v- v J ini7 P very short variability time as due interstellar scm- 



peratures in the radio regime which may exceed lO'^^ K, or 
requires relativistic beaming with extremely high Doppler 



tillation requiring the angular diameter of the most compact 

component to be ~ 6 micro-arcsec, or smaller, correspond- 

ractors, coherent radiation mechanisms, or special geometric . , ,. , , i j_ i i i rr^-ir r-i 

„ /-,-,, „ , „ n , . , n . , mg to a solid angle subtended by the IDV core il — ih 

ettects (Wagner & Witzel 1995). buch high brightness tem- , „ „ „„ ,r,-22 mi r i- r ii j_ j_ i n 

^ ° , , , „ , , , „ where S2i ^ 6.79 x 10 sr. The traction ot the total flux 

peratures would be well above the Compton catastrophe ,, , , , , . . i -.i . , rr^^r 

, 11 ^ they estimated to be associated with this compact IDV core 

limit Tb < 10 K imposed by mverse-Compton scattering .„ i-i-w 

, ^. . , , IS 6c ~ 0.15, and the corresponding brightness temperature 

(Kellermann & Paulim-Toth 1969, blysh 1992, Kardashev ^ . „ „,t ■ m r 1^^l4/r-. /<-.\ a i^i i i i 

^ , , , , , , at 4.8 GHz IS Ts « 5 X 10 (rii/ri) K. Although much lower 

2000) when the photon energy density m the emission region „ , .i . , . i . . . .,i ,i 

' n 1 1 n 1 than J var, to reconcile this brightness temperature with the 

exceeds the energy density m the magnetic held. See Krich- i^iit^i- i i ^ ■ i i 

, , s , t/ n N r ~ 10 K limit would appear to require a very large Doppler 

baum et al. (2002) and Kedziora-Chudczer et al. (2001) tor r , n -,^3 i i i. . m 

^ tactor, V ~ 10 , or an even larger angular diameter, lo 
recent reviews ot IDV sources. , i ^ ^ 1^3 ^^ 

obtain a Doppler factor as large as 10 would require not 
The radio-loud quasar PKS 0405-385 is an extreme ex- , ■ , r r ^ r-n^r-./r, mr, 1 ^ 1 

only a jet Lorentz tactor 01 1 > D 1 = 500, but also very 
ample or an intra-day variable source with variations on , ,. ^ r j_i ■ ^ • ^ i- r • 1 j_ / -^.i • 

. , ^ r, -, n /T^ , . ^, , ,1 „^„N close alignment 01 the jet axis to our line ot sight (withm 

timescales ot fmv ~ 0.1 d (Kedziora-Chudczer et al. 1997). 1 /n , n ion mi 1 -i-i r 1 i- 

, ; ~ 1/1 < 0.1 ). i he probability ot such an alignment occur- 

Makmg the assumption that the emission region subtends . , , . 1 ^-r•2 ^ 10-6 • j_i ■ 1 ■ 1 

, tnt,^ n\2/j2 u n • T-> 1 ring by chance IS then ~ l/4r < 10 , m this case, which 

solid angle ~ -Ki^^.octiDV L>) /"«, where U is the Doppler 1 1 ■ 1 1 r ^ -i -i-^ 

" \ / ; w ff makes the very high Doppler factor possibility unattractive, 

factor and dg is the diameter distance to the source, one can 

The effective distance to the interstellar scintillation 
screen is crucial in determining the angular size of the 

* email: rprother@physics.adelaide.edu.au source, and hence its brightness temperature. The dis- 
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tance used by Kedziora-Chudczer et al. (1997) and Walker 
(1998) is effectively the scale-height above the Galac- 
tic plane of electron number density squared, i.e. Z2 = 
zn1{z)dz/ n^{z)dz, for which Z2 ~ 500 pc in the 
model of Taylor and Cordes (1993) for the free electron 
distribution in the Galaxy. However, that model was de- 
veloped mainly for the consistent determination of pulsar 
distances from dispersion measures, and is most accurate at 
low galactic latitudes where the majority of pulsars are ob- 
served. Taylor and Cordes (1993) themselves warn that one 
should be aware of uncertainties in their model associated 
with this. The column density, Uedz, is fairly accurately 
determined, but the scale height of electron number density, 
zi = zne{z)dz/ ne(z)dz, is less well determined, and 
Z2 is even less accurately known. For example, in the recent 
model of Gomez et al. (2001) for the free electron distribu- 
tion in the Galaxy Z2 ~ 300 pc. It is interesting to note that 
Beckert et al. (2001) suggest typical distances to the scat- 
tering medium of 200 pc, and that a scale height of about 
100 pc seems to be required to explain IDV in the case of 
0917-^624. 

Krichbaum et al. (2002) mention the possibility of an 
extremely clumpy ISM. If this is the case, one could well 
question the use of an effective screen distance as it could 
be that the scintillation is due to an individual ionized cloud 
much nearer to us than the average distance. In fact, if the 
distribution is highly peaked in the Galactic plane the most 
probable, rather than average, distance is small. The dis- 
tance to the ionized cloud responsible for the extreme scat- 
tering event in 0954-1-658 is estimated to be ~ 150 pc (Cimo 
2002), while Dennett- Thorpe & de Bruyn (2000, 2002) sug- 
gest that the scattering region for IDV in J1819+3845 may 
be located at about 20 pc, and possibly associated with 
the local bubble. Hence, I believe there could be consider- 
able uncertainty in the distance to the scintillating material 
responsible for the extreme IDV in PKS 0405-385. If, for 
example, the distance were smaller by a factor of 5, this 
would translate into a factor of 25 reduction in the bright- 
ness temperature. A brightness temperature of ~ 2 x 10^^ K 
at ~ 5 GHz for the IDV core in PKS 0405-385 is still high, 
but I will show that it can be achieved by standard electron 
synchrotron radiation using quite moderate Doppler factors 
if one takes into account possible geometries of the emission 
region. 

In a recent paper (Protheroe 2002) I have explored the 
effect of emission region geometry on flux variability, and 
on the relationship between observed intensity and energy 
density for various source geometries for the case of optically 
thin emission and found that the average energy density in 
the source can be much less than one would estimate simply 
from the observed intensity. Although the radio emission 
from IDV sources is usually assumed to be optically thick, 
if this is not the case then the above result may also have 
important implications for IDV sources as the photon energy 
density responsible for causing the brightness temperature 
limit may actually be a few orders of magnitude lower than 
estimated from the intensity. In that case, lower Doppler 
factors would be required to avoid the Compton catastrophe. 
In this paper I shall explore the parameter space, including 
emission region geometry, of models able to reproduce the 
observed radio emission of the IDV core of PKS 0405-385, 
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Figure 1. Observed fluxes of PKS 0405-385 divided by are 
plotted as gray symbols and vertical lines. Contemporaneous VLA 
and ATCA data, multiplied by Sc, are also re-plotted as black 
symbols and represent the flux from the IDV core. The right- 
hand axis shows the brightness temperature assuming the IDV 
core subtends solid angle Oi . Data; pluses — non-contemporaneous 
data obtained from the NED archive; triangles and squares - 1996 
and 1998 VLA/ATCA data (D. Jauncey personal communication 
1999); vertical solid lines - 1993-4 SEST data (Tornikoski et al. 
1996); dashed vertical lines at 230 GHz and 4.3 X 10^"' Hz - (S. 
Wagner personal communication 1999). The solid curves are fits 
for synchrotron radiation by mono-energetic electrons with vi = 
lO^'i Hz and Uc = lO" Hz (left) or Uc = 10^^ '^^ Hz (right). 

and I shall model its spectral energy distribution (SED) from 
radio to gamma-ray frequencies. 



2 FITTING THE SPECTRUM OF PKS 0405-385 

In Fig. I plot the available flux measurements of 
PKS 0405-385, divided by frequency squared, as gray sym- 
bols and gray vertical lines. Contemporaneous VLA and 
ATCA data (squares and triangles in the 1.4 GHz-43 GHz 
range) are also re-plotted as black symbols at the observed 
fluxes multiplied by Sc, the fraction of the total flux assumed 
by Kedziora-Chudczer et al. (1997) to be associated with the 
IDV core. The right axis shows the brightness temperature 
of the IDV core inferred by Kedziora-Chudczer et al. (1997) 
assuming the IDV core subtends solid angle fii. It is only 
the IDV core component of these data that I am concerned 
with fitting. Nevertheless, it is interesting for comparison 
to include non-contemporaneous data at other wavelengths 
whose origin may or may not be associated with the same 
emission region. 

Over the 1.4 GHz-43 GHz range, the intensity can be 
fitted well by oc i/^''^ (Tb oc i/"^/^). Such a spectrum 
would occur naturally if the emission were optically thin and 
if this frequency range were well below the critical frequency 
of the lowest energy electron. Note that for an XSelectron 
with Lorentz factor 7 the critical frequency is 

i^c = ^J-^ = 4.19 X 10%' Bx Hz (1) 
AnmeC 

where e is the electron charge (statcoulombs), B± is the com- 



The brightness temperature problem in extreme ID V quasars: a model for PKS 0405-385 3 



ponent of magnetic field (gauss) perpendicular to the elec- 
tron velocity, and rrie is the electron mass (grams). Such an 
electron spectrum could occur naturally if the electrons were 
produced as secondaries of other particles, e.g. Bethe-Heitler 
pair-production by relativistic protons, or if explosive re- 
connection were responsible for electron acceleration. For 
the purposes of this paper, I shall adopt a mono-energetic 
electron distribution. 

For a mono-energetic isotropic electron distribution the 
synchrotron emission coefficient is given by 



47r 



-F{x) 



(2) 

where Tie is the electron number density (cm~'^), x = v/pc 
and F{x) = x ^5/3(C)'^C (s^e e.g. Rybicki and Lightman 
1979). At low frequencies v ^ i/c 



F{x)^ 



(f) 



(f) 



1/3 



(3) 



The absorption coefficient for an isotropic electron dis- 
tribution N{E) is 



dE P{iy,E)E^ 



N{E - hu) N{E) 



(4) 



For the mono-energetic electron distribution considered here 
N{E) = neS{E — 7mc^), and assuming hv <^ E, I find 



d 



33/25x7^ dx 



F{x) 



(5) 



and this is plotted in Fig. |21 Note that at low frequencies 



32n erie 



27Bx7^r 



181.4neBx''S 



21/3 

/3, 



-5/3 



5/3 



(6) 

and this is plotted as the dotted line in Fig. |5| From equa- 
tions JHJ and Q the source function at low frequencies is 

3 2 



(7) 



which is identical to the Rayleigh- Jeans approximation to a 
black body spectrum of temperature 



T = 



3'ymei 



(8) 



Thus, in the optically-thick very-low frequency range the 
brightness temperature is constant at Tb ~ 4.45 x 10^7 K, 
and then there is a transition to Tb oc v'^^'^ at frequency ui 
where the optical depth Tv = J a^di is unity, i.e. r^(i'i) = 1. 
Hence, TBiyx) can be used to estimate the Lorentz factor of 
the electrons. 



As noted earlier, the data appear to show 1^ oc 



,1/3 

over the 1.4 GIfz-43 GHz range corresponding to the con- 
temporaneous IDV observations. At other frequencies it is 
not known whether the observed emission is due to the same 
compact component, or is due to a larger region, perhaps 
farther along the jet. However, data at other frequencies 
may still be used to constrain the models. Since the bright- 
ness temperature at low frequencies is proportional to the 
electron Lorentz factor, the brightness temperature problem 
is minimized by using the lowest possible Lorentz factor. I 
therefore take the highest frequency v\ which is just con- 
sistent with the 1.4 GHz data, and adopt vi = 10^'^ Hz. 




Figure 2. The synchrotron absorption coefficient vs. x = v/vc 
for mono-energetic electrons with Lorentz factor 7 and number 
density ric (cm"^) in a magnetic field with perpendicular compo- 
nent _Bx (gauss). The dotted line shows the low- frequency limit 
ss 181.4neB^''^7-'5/3i/-5/3^ 



Apart from the normalization, for which 1 choose to take 
Fx = F^{vx) ~ 8.9 X 10"^^ erg cm^^s^^Hz"^ (after multi- 
plying by Sc ~ 0.15), the other parameter determining the 
fit is the value of the critical frequency Vc . The lowest critical 
frequency Uc which is consistent with the 43 GHz radio data 
is Vc « 10^ Hz, and the highest critical frequency consistent 
with the R-band data is Vc ~ iq13-65 jj^^ ^-^^ j shall adopt 
these two values in my modeling. The resulting brightness 
temperature fits are shown by the solid curves in Fig.Q The 
challenge is now to find a combination of source parameters, 
i.e. Doppler factor, jet-frame magnetic field, electron num- 
ber density and emission region geometry, that gives the 
required intensity, and is physically possible. 

If one views the emission region along the jet axis, 
one obtains the largest Doppler boosting for any given jet 
Lorentz factor, and I shall assume this to be the case for 
PKS 0405-385. This may well be justified as such extreme 
IDV radio quasars are extremely rare. Generally, the lower 
the energy density of synchrotron radiation photons in the 
emission region, the lower will be the Doppler factor re- 
quired to avoid the Compton catastrophe on synchrotron 
photons. However, as I shall show, very large Doppler fac- 
tors can cause a Compton catastrophe on cosmic microwave 
background radiation (CMBR) photons. I shall discuss the 
dependence of synchrotron photon energy density on the ge- 
ometry of the emission region in the next section. 



3 SYNCHROTRON PHOTON ENERGY 

DENSITY AND AND EMISSION REGION 
GEOMETRY 

In a separate paper (Protheroe 2002) I have discussed the 
influence of emission region geometry on photon energy den- 
sity in the emission region for the optically thin case. Here 
I shall extend this work to include optical depth effects and 
use a cylindrical emission region. I shall consider the case 
of a cylinder of length I and radius r, with uniform emis- 
sivity ju and absorption coefficient a^, and determine the 



4 R. J. Protheroe 



photon energy density (U,^) averaged over the volume of the 
cy Under. 

At some arbitrary point r within the cyhnder the inten- 
sity from direction defined by unit vector e is 



Ii,{f, e) = —{1 - exp[~a^x{f, e)]} 



(9) 



where x{r, e) is the distance from f to the boundary of the 
cylinder in direction e. The energy density at r is 



(10) 



Using a Monte Carlo method one can sample a large number, 
Nd, of directions ei,i = 1, . . . ,Nd, distributed isotropically, 
and then set 



(11) 



The energy density averaged over the volume V — nr^l 
of the cylinder is then 



(12) 



Using a Monte Carlo method one can sample a large num- 
ber, Np, of points ffcjfc — 1, . . . ,Np, distributed uniformly 
throughout the volume of the cylinder, and then set 



Hence, 



NdNp c 



(13) 



- exp[-a^x{fk,e,)]} , (14) 



and in this way one can calculate (Uv) for various £/r values, 
and absorption coefficients a,/. Of course for emission in the 
jet frame, all the variables in equation 1141 would be jet- 
frame variables. 

While the average energy density inside the cylinder is 
fixed for any set of j^,av,£/r, and r, the intensity observed 
when viewing the cylinder can depend strongly on viewing 
direction. In order to obtain the highest observed intensity 
one would look in a direction such that the projected area 
of the cylinder is smallest. Making the axis of the cylinder 
coincident with the jet axis, this would be achieved if one 
viewed the emission region at 9' = n/2 if the cylinder was 
short {i' /r <C 1), and for ^' = if the cylinder was long 
i£'/r 3> 1), where primed coordinates correspond to jet- 
frame variables. I shall consider the latter case, i.e. viewing 
the emission down the jet axis {9' — 0). The jet- frame inten- 
sity is then simply given by II, = {jl,/a'^,)[l - exp(-T^,)] 
where r^; = a'^ii' . 

In Fig. 01 plot the ratio {Ul) /{A-kII/c) against r^, 
for various values of I' /r. As can be seen, the effect is 
very important where the emission is optically thin, e.g. for 
I'/r = 10^ the average energy density is almost a factor 10 
lower than would be expected from the observed intensity 
(note that the optical depth is Lorentz invariant r^/ = t^). In 
the next section I shall consider various i'/r values when ex- 
ploring the parameter space which could apply to the emis- 
sion region for IDV in PKS 0405-385. 




Figure 3. Average energy density inside the emission region di- 
vided by Att/c times the observed intensity vs. optical depth for 
various £'/r. 



4 MODEL PARAMETERS 

The fits shown by the solid curves in Fig. correspond to 



Fi = F^{vi) « 8.9 X 10"^^ ergcm-^s-^Hz-\ 



W "^ Hz, 



and Vc = 10^^ Hz or IQiS-Ss Hz. In this section I shaU deter- 
mine the combinations of physical parameters which could 
give rise to these values, assuming a cylindrical emission re- 
gion geometry, an isotropic mono-energetic electron distri- 
bution in the jet-frame, and observation along the jet axis. 
The physical parameters describing the emission region are 
the perpendicular component of magnetic field B'^ as mea- 
sured in the jet-frame, the Doppler factor D, the ratio I' /r of 
the emission region cylinder length to its radius as measured 
in the jet-frame, and an equipartition factor which gives the 
ratio of relativistic electron energy density to magnetic en- 
ergy density 



7? = (7 n^nieC ) [ - — 



(15) 



where the (3/2) factor arises because B'^ is the perpendicu- 
lar component of magnetic field, rather than its magnitude, 
and 7' is the jet-frame electron Lorentz factor. 

The solid angle subtended by the source is obtained 
from the angular radius of the source, i.e. Q. — 1^91 f^nd the 
radius of the cylindrical emission region is obtained from the 
diameter distance and angular radius, i.e r — dg9i^2- Note 
that for PKS 0405-385 at redshift z = 1.285 the diameter 
distance is dg = 5.68 x lO'^^ cm for a ACDM model with 
Qru = 0.2, = 0.8 and Hq = 70 km s"^ Mpc"^ (here 
Q specifies the fraction of the cosmological closure density, 
whilst elsewhere in this paper is used for solid angles). 

Noting that Ii, = F^/fl, and that in the optically thick 
region of the spectrum I^, — j^/a^, from equation Q one 
may obtain the Lorentz factor the electrons would have if 
the emission took place in the observer frame, 71, and hence 
find the jet frame Lorentz factor 



(16) 



Next, using equation Q and that, by definition, the optical 
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depth at jet-frame frequency v'l — i/i/D must be unity 

181.4 (//r)(6li/2d9)niB'l''^ x 

/ 2Fi V''^'^^ 
[sD^el^^m^uf ) [d) 

Substituting for nj. from equation l|15^ with 7' from equa- 
tion 1II6I I and solving for B'j^ I obtain 

d' a no in23r-i-13/8 -11/8^-19/8 „ 

B± — 4.98 X 10 D ' B^i^ X 

xFid-'/* G. (18) 

Then, using equation Q for the jet-frame critical frequency, 
I'c = ^cjD, I obtain 

Bi = 4.38 X 10-'*°D!^iV,6»^/2-F'r^ G. (19) 

Finally, solving simultaneous equations I18II and 1191 I ob- 
tain 

D = 4.37xlO^^,-'^'Vr'^''2^ X 

x^-///V^X-/^(^4)"'' (20) 

Bl = 1.91xlO-^Vf/^Vf/'^<)f X 

xf-/^-/^ ' G. (21) 

Fig. m shows the magnetic field-Doppler factor param- 
eter space for models fitting the radio intensity of the IDV 
core of PKS 0405-385. Equations and define mod- 
els which will give synchrotron spectra of mono-energetic 
electrons determined by F\, v\ and Vc, and these models are 
represented as solid lines in Fig. |1] corresponding to either 
a fixed value of {rjl' /r) while varying 9i/2 as a parameter 
in these equations, or a fixed value of 61/2 while varying 
(t\(! jr). All other variables are fixed at the values appropri- 
ate to PKS 0405-385, with Fig. El; a) being for = 10" Hz 
and Fig. Hb) being for = 10^=* '^^ Hz (yx = 10^'^ Hz in 
both cases). 

If one allows the relativistic particle energy density to 
exceed the magnetic energy density, then very low Doppler 
factors are possible. However, if the equipartition factor is 
greater than unity, this would be unstable. It is therefore 
probably unrealistic for the relativistic particles to be too 
far from equipartition with the magnetic field, and so one 
should perhaps take more seriously results corresponding to 
rj — 1, which I shall assume in what follows. Nevertheless, 
IDV is truly a time-dependent problem, and episodes with 
r? > 1 are not ruled out completely. 

4.1 Avoiding Compton catastrophes 

If electrons are injected with jet-frame Lorentz factor 7' then 
the presence of dense radiation fields provides target pho- 
tons for inverse Compton scattering, and the generation of 
components in the SED at X-ray and gamma-ray frequen- 
cies which may or may not exceed observed X-ray fiux and 
gamma-ray limit. If the electrons are accelerated in a quasi- 
continuous process such as diffusive shock acceleration by 
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Figure 4. Magnetic field— Doppler factor parameter space for 
models fitting the radio intensity of the IDV core of PKS 0405-385 
with 1/1 = 10^1 Hz and (a) Uc = lO" Hz or (b) Uc = 10i3-65 jj^, 
i.e. giving the left and right solid curves in Fig. 1. The two sets 
of solid curves labeled by log(»?^'/^) ^^d log(6i/2) ^re given by 
equations I20i and I21i plotted using 61/2 and {r]£' /r), respec- 
tively, as the parameter, and are used to obtain the values of 61/2 
and {rji' /r) for any combination of B± and D that are necessary 
to give the required radio spectrum. The shaded region on the 
right shows where the energy losses are dominated by IC on the 
CMBR and corresponds to equation 1291 . and the shaded regions 
on the left show where the energy losses are dominated by IC on 
synchrotron photons and correspond to equation 1301 for various 
values of log{£' /r) indicated. Fifteen models, i.e. combinations of 
B±, D, Q\i2 '^nd ijil' /r), which will be discussed later are labeled 
A to O. 



non-relativistic or mildly-relativistic shocks, then the pres- 
ence of dense radiation fields may lead to excessive energy 
losses and prevent the electrons reaching the Lorentz factor 
required to fit the observed radio spectrum. In either case, 
these effects may occur wherever the photon energy density 
in the emission region becomes comparable to or exceeds the 
energy density in the magnetic field. 
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In addition to the synchrotron emission itself, I shall 
assume that emission region is sufficiently far along the jet 
that the only other non-negligible field for inverse-Compton 
scattering is the cosmic microwave background radiation 
(CMBR) . If the emission region were close to the central en- 
gine then accretion disk radiation, broad line cloud emission, 
and torus infrared emission could also play a role (Donea & 
Protheroe 2002). 

The average jet-frame energy density of the synchrotron 
emission inside the emission region can be estimated for 
given i' /r values from the observed intensity correspond- 
ing to v^F,^ plotted as the solid curves in Fig.Q The energy 
density in the jet frame is given by 



where 



4tv 



— / M[£'/r,T4i^)]L{,,,er/2)diy, 



(22) 



(23) 



r^(i/) ^ (fi/i^f^^ from Fig. El and M[£' /r,T^{iy)] = 
(i7^/)/(47r/^//c) given in Fig. |3] The main contribution to 
the energy density for the spectrum fitted to the radio ob- 
servations of the IDV core in PKS 0405-385 comes from fre- 
quencies near to Uc since this is where the SED, i/Fv, peaks. 
This spectrum is optically thin near Vc and so I may make 
the approximation 



47r 



M{e'/r,T, < 1)— / 01/2 )di/ 



M(//r,T. « l)ff(/.i,i/e)i^ie-/2, 



where 



du. 



(24) 



(25) 



F^{vi) 

For the spectra shown by the solid curves in Fig. 
M(l,r^ < l)g(10^-^Hz, 10"Hz) = 45.14 cm~^sr Hz and 
M(l,Ti.<l)ff(10''-iHz,10^^ ''5Hz) = 1.47 x lO'^ cm-^sr Hz. 

In the rest frame of the host galaxy the CMBR would 
have had temperature (1 -I- z)To, where To = 2.735 K, and 
would have been isotropic. I shall distinguish between the 
Doppler factor for boosting from the jet frame to the host 
galaxy frame, 8, and the Doppler factor for boosting from 
the jet frame to the observer frame, D — 5/(1 -I- z). In the 
jet frame, the black body temperature would depend on the 
angle ip' with respect to the jet axis at which an observer in 
the jet frame looks 



T'(^') = 



(1 + z)To 
r(l - /3j cost/)') 



(26) 



where f3jC is the jet velocity as measured in the host galaxy 
frame, F = (1 — /3^)~^^'^ is its Lorentz factor, and is given by 
r = {5^ + l)/(2(5) for the case of observation of the AGN jet 
at angle 9 — to its axis. The energy density of the CMBR 
in the jet frame is then 



CMBR — 



(1 + zfaT^f f dn' 

47r / [r(l-/3jcos6i')l* 

^ It 



z)'^aTo. 



(27) 



where a = 4(t/c and a is the Stefan-Boltzmann constant. 
Note that whereas the jet frame energy density of syn- 
chrotron radiation decreases with assumed Doppler factor, 
the energy density of the CMBR increases with assumed 
Doppler factor. This means that the Compton catastrophe 
can not be avoided by using an arbitrarily large Doppler 
factor. In fact, very high Doppler factors will result in a 
Compton catastrophe due to inverse-Compton scattering on 
the CMBR. This has important consequences which I shall 
address next, and in subsequent sections. 

For efficient synchrotron radiation to take place the en- 
ergy density in the magnetic field should exceed the energy 
density in the radiation field 



3 B'i 



> 



D-^{U.y4£'/r)) + -r\l + z)\T,t 



(28) 



Since the synchrotron photon energy density is large in the 
regime of low Doppler factors, while the CMBR energy den- 
sity is large in the regime of high Doppler factors, I shall 
obtain separately the minimum magnetic field required in 
each case. For PKS 0405-385 being viewed along the jet 
axis, and using the approximation F ~ (1 -|- z)D/2, valid for 
F ^ 1, the jet-frame energy density of the magnetic field is 
less than that of the CMBR when 



< (4/3)^/^(l + zfriD 



(29) 



and this is independent of all model parameters except D. 
This case is shown by the shaded area on the right in 
Figs.Ha)-(b). 

In the case of synchrotron photons as targets for in- 
verse Compton scattering, i.e. the synchrotron self-Compton 
(SSC) process, inverse Compton losses dominate when 



< 



4(^/3)^/^D-^{f/.,„(//r,ei/2))i/^ 



(30) 



Because in this case the photon energy density depends on 
9i/2 and I' /r, in addition to D, the boundary between mod- 
els {D,B'^) is parametrized by substituting D from equa- 
tion (I19II into equation 1301 1 and solving for B'^ as a function 
of the parameters, and then using equation (I19II to obtain 
_D as a function of the parameters. The resulting minimum 
value of B'j^ is plotted against D in Figs.2Ia)-(b) for various 
values of £' /r as the dotted lines (labeled by £' /r) bounding 
the shaded areas on the left. 

The gyroradius, r'g = "f' c/ [B'^lub) must be much 
smaller than the radius r of the cylinder and so this pro- 
vides, in principle, an additional lower limit to B'^. However, 
this limit is well below the lower limit from equation 1281 
already plotted, and therefore does not further constrain the 
models. 

Fifteen potential models which will be discussed later, 
i.e. combinations of B^, D, 6^/2 and (r\£! jr), are labeled 
as A to O in Fig. |1| Before proceeding, we should check if 
any of these models is ruled out by being optically thick to 
Thomson scattering. The Thomson optical depth is tt = 
£!ri^aT where ctt is the Thomson cross section. Taking 
from equation 11511 with 7' from equation 1161 I obtain 



TT = 2.11 



X V<S^^B'\De{ 



/2 



1)- 

r 



(31) 



For the fifteen potential models the Thomson optical depth 
ranges from tt = 7.04x 10"^^ (Model A) to tt = 1.83x10"'^ 
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(Model O), and so Thomson scattering can be neglected in 
this case. 



where 

<Wsyn(l^0)> = T-T V (''O,^' = 0) 



(38) 



5 CALCULATING THE SPECTRAL ENERGY 
DISTRIBUTION 

I shall consider the following three emission processes: syn- 
chrotron radiation (syn), inverse Compton scattering of syn- 
chrotron photons (SSC), and inverse Compton scattering of 
CMBR photons (ICM) . Assuming the inverse Compton scat- 
tering is in the Thomson regime, the rate of energy loss by 
either process is 



dE' 



(32) 



where U' = U'g for synchrotron losses, U' = (U'syyi) for SSC 
losses and U' = f/cMBR for inverse-Compton losses on the 
CMBR. The fraction of emitted radiation by each process is 
then 



/syn 

/ssc 

/iCM 



(JJly^)/{U'B + (f/s'yn) + t/cMBR.) 
UcmBkKU'b + (JJ'sy-n) + t^CMBft) 



(33) 
(34) 
(35) 



This does not mean, however, that the observed energy flux 
for the inverse-Compton scattered CMBR is in this ratio to 
the other two components because, whereas the synchrotron 
emission and SSC emission is isotropic in the jet frame, the 
inverse-Compton scattering of the CMBR is not because of 
the anisotropy of the CMBR in this frame. One may well 
expect strong peaks in the SED due to inverse-Compton 
scattering for models close to or within the shaded areas 
in Figs. |3Ia)-(b). Slysh (1992) also noted that gamma-ray 
emission may occur in IDV sources. 



5.1 Synchrotron radiation 

The intensity of synchrotron emission in direction ^ = in 
the observer frame is simply given by 

/r (J^, e = Q)= D^Sl, {u/D){l - exphr^, {,./D)]} (36) 
where tI,[v/D) = ct'^,{v/D)l'. 



5.2 Inverse Compton scattering of synchrotron 
photons (SSC) 

In the present paper I shall assume that we view the emission 
emitted at angle 6' = 6 = Q to the jet axis. The electrons 
are assumed to be isotropic and mono-energetic in the jet 
frame, and have Lorentz factor 7' = 11/ D, and velocity 
P'c = (1 — l/7'^)^/^c. For simplicity, I shall also take the 
synchrotron target photons to be isotropic in the jet frame, 
while using the results from Section 3 to normalize their 
spectrum. 

For isotropic target photons in the jet frame with fre- 
quency Uq, the frequency of the scattered photons range be- 
tween and 47''^//o. The jet-frame SSC emissivity (erg cm~'^ 
s~^ sr~^ Hz~^) is then 



./SSC n'^ihv'cjTC 



47'- 



'^KyM))hc{x) (37) 



is the average jet-frame synchrotron specific photon num- 
ber density (photons cm"'' sr~^ Hz~^), x — v' l^i v'o and 
fic{x) = 22;lnx + 3:: + l-2a::^ (Blumenthal and Gould 1970), 
giving 



tI SSC, I qI n\ 1)1 -'SSC 



(39) 



Finally, this is Doppler boosted to the observer frame 
l^^'^{u-e = 0) = D^l'Jf^^{v/D;e' ^Q) (40) 



5.3 Inverse-Compton scattering of the CMBR: 
gamma ray production 

In the jet frame, the black body temperature would depend 
on the angle 6q with respect to the jet axis at which the 
CMBR target photons propagate 



T'{e',) = 



(1 + z)To 
V{1 + PjCosO'q) 



(41) 



such that their specific photon number density (photons 
cm^"' sr^^ Hz^^) is 

ncMBR(i'o,^'o) = — J7— rTTTFTT^TTi 7- (42) 



exp[KAr'(eo)]-i' 

For target photons in the jet frame with frequency v'q 
propagating at angle Q'^ to the jet axis, the frequencies of the 
photons scattered by electrons propagating parallel to the jet 
axis are uniformly distributed between and i'max(i'oi ^0) = 
27'^i/o(l — cos Sg) in the approximation that the scattered 
photons are isotropic in the electron rest frame. The jet- 
frame emissivity (erg cm~^ s~^ sr~^ Hz~^) for IC on the 
CMBR in the jet direction is then 



./ICM;-/)/ n\ 

3v> (d =0) 



/, / 1"^ dcose'Q I , 

n^niy arc j (1 — p cosc^q) 



, / ' 

dvo- 



where v''o^^"^{Oq) ■ 

:0) 



11"=^ {v; 9 = 0) 



v' I2i\\ - /3' cos 6*0). Finally, 

= ^X!^^(e' = o), 

,3 r'lCM^^/^.g/^Q) 



(43) 



(44) 
(45) 



6 DISCUSSION 

The angular radius inferred by Kedziora-Chudczer et al. 
(1997) is 6*1/2 « 1.5 X 10"" rad, i.e. log(6ii/2) = -10.83. 
Examining the parameter space for Vc = 10^' in Fig. 21 a), 
one sees that all models with this angular radius lie well in- 
side the region where IC on the CMBR dominates the energy 
losses of electrons. For r\ = \ (equal energy density in mag- 
netic field and relativistic electrons) and £'/r = 1 this would 
correspond to a point roughly mid-way between models A 
and B, and would require a Doppler factor of D ~ 5000. By 
choosing jr > 1 one may reduce the Doppler factor, but 
unrealistically large (! jr values are needed to reduce D by 
a large factor. To reduce the Doppler factor to a reasonable 
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Figure 5. SED of IDV core of PKS 0405-385 obtained from ob- 
served fluxes assuming flux fraction Sc only for the ACTA/VLA 
data as in Fig. 1. Data additional to those shown in Fig. 1: X-ray 
data at ~ 10^® Hz - 1997 ASCA fluxes (L. Kedziora-Chudczer 
personal communication 1999); gamma-ray data — EGRET up- 
per limit (Hartman et al. 1999). Sensitivities of the IBIS instru- 
ment on INTEGRAL (Si), GLAST (Sq) and HESS (Sh) indi- 
cated (the expected sensitivity of CANGAROO III is similar to 
that of HESS). The calculated SED is shown for various models 
with the synchrotron component as the solid curve on the left, the 
SSC component in the middle (when present) and the scattered 
CMBR component on the right, (a) The SED for models A-E 
in Fig. Ilia) with Uc = lO^ ^, rj = 1 and {£' /r) = 1: A - dotted 
curves, B — short dashed curves, C - chain curves, D — triple-dot 
dashed curves, E - long dashed curves. 

value, i.e. D < 100, would require the angular diameter to 
be larger by only a factor of ~ 5 which, as discussed earlier, 
I believe is quite possible (for an angular diameter larger by 
a factor of 10, even D ~ 10 is possible). 

6.1 Spectral energy distribution 

Because of the dominance of IC losses on the CMBR for 
models A-D, one would expect the SSC peak at gamma- 
ray energies in the SED to exceed the synchrotron peak for 
these models. The SEDs have been calculated for models A- 
E and are shown in Fig. IHa) and, as expected, show strong 
gamma-ray emission for models A-D. Non-contemporaneous 
X-ray and gamma-ray data are shown, and may be used with 
caution as upper limits. I have added to Fig.|^the sensitiv- 
ity of the IBIS gamma-ray detector on INTEGRAL (Par- 
mar et al. 2002), and the expected sensitivities of GLAST 
(Gehrels and Michelson 1999), and HESS (Hofmann et al. 
2001) and CANGAROO III (Enomoto et al. 2002), the last 
two being southern hemisphere atmospheric Cherenkov tele- 
scopes nearing completion. Taking the X-ray data as an up- 
per limit does not rule out any of these models. However, the 
EGRET upper limit obtained from figure 3 of Hartman et 
al. (1999) appears to exclude models A-C but be consistent 
with models D-E (note, however, that the EGRET data are 
not contemporaneous with the IDV data). Models D-E have 

^ 6 X 10""-10"^° rad and D ^ 180-40, respectively. 

As expected from equation 1201 increasing Vc from 



-5 't' / n 

16 t 
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log(i//Hz) 
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16 t 




log(i./Hz) 

Figure 5 - continued (b) The SED for models F-J in Fig. IHb) 
with i/c = 10^^-^^, ri = 1 and {£' /r) = 1: F - dotted curves, G 

— short dashed curves, H — chain curves, I - triple-dot dashed 
curves, J - long dashed curves, (c) The SED for models K-O in 
Fig.EJb) with = 10l3-65^ r; = 1 and {I' /r) = 316: K - dotted 
curves, L — short dashed curves, M - chain curves, N - triple-dot 
dashed curves, O - long dashed curves. 

10" Hz to 10" 65 Hz reduces D by a factor of 10 if one keeps 
all other parameters fixed (note, however, that this reduction 
in D is accompanied by an increase in B'j_). This is clearly 
seen in Fig. Elb) which is for Uc = IQ^'^-^^ Hz, and thereby 
presents more opportunity for obtaining reasonable values of 
D. For r] = 1 and £' /r = 1, models F-J have the same angu- 
lar radius range (10""-10"^° rad) as models A-E, but with 
considerably lower Doppler factors, and these models range 
from being in the IC on CMBR dominated regime (model F 

- expect highest gamma-ray fiux, lowest SSC fiux) to being 
in the SSC dominated regime (model J - expect highest SSC 
fiux, lowest gamma- ray fiux) . The calculated SEDs for mod- 
els F-J are shown in Fig.|3^b) and show the trends expected. 
Possibly models I and J may be ruled out by the X-ray data 
taken as an upper limit (note, however, that the X-ray data 
is not contemporaneous with the IDV data). Models F-I 
predict high gamma-ray fiuxes from MeV energies to 100 
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GeV energies, depending on model, which could be tested 
observationally by INTEGRAL, GLAST, HESS and CAN- 



'1/2 



10" 



-6 X 10" 



rad 



GAROO III. Models F-I have 
and D ~ 1000-20, respectively. 

As also expected from eauation l20l increasing £' /r from 
1 to 316 results in a reduction in D by a factor of ~ 2 
(plus a reduction in B'^) if one keeps all other variables 
fixed. Models K-O in Fig. IHb) have all other parameters 
the same as models F-J, and the resulting SEDs are shown 
in Fig.l^c). Because of the reduction in synchrotron photon 
energy density associated with £' /r ^ 1, although not iden- 
tical, the SEDs of models K-O are qualitatively similar to 
those of models F-J, with model O being ruled out by both 
the gamma-ray limit and the X-ray data, although again I 
note that these data are not contemporaneuos with the IDV 
data. Models K-N have ^i/a ~ 10""-6 x 10"" rad and 
D « 400-8. Model M, having 6*1/2 only ~ 5 times that as- 
sumed by Kedziora-Chudczer et al. (1997) has a Doppler fac- 
tor of ~ 30, and lower Doppler factors are possible with, for 
example even larger £' /r. Bearing in mind that no attempt 
has been made to fit the (non-contemporaneous) X-ray and 
gamma-ray data, and that a simple mono-energetic electron 
spectrum has been used, it seems that models which are able 
to fit the IDV radio data tend to predict observable fiuxes of 
X-ray, and/or gamma ray emission at sub-GeV and/or ~10 
GeV energies. Use of a more sophisticated electron spec- 
trum is unlikely to alter this conclusion. Hence, it may well 
be profitable for space gamma ray telescopes such as INTE- 
GRAL and GLAST, and southern hemisphere atmospheric 
Cherenkov telescopes such as HESS and CANGAROO HI 
to look for emission from PKS 0405-385 and other extreme 
IDV radio galaxies. 



6.2 Variability 

The rapid intra-day variability in PKS 0405-385 is almost 
certainly due to interstellar scintillation effects. The radio 
fiux does, however, also vary on longer timescales - 10 per 
cent change in 2 months in July 1996 (corresponding to ~ 
1 yr for 50 per cent change), and has very intense for periods 
of about a month (Kedziora-Chudczer et al. 1997). 

A mechanism that may cause such variability would be 
the emission region moving along the jet and passing a region 
where external factors cause compression/expansion of the 
emission region, increasing/decreasing the magnetic field, 
and causing adiabatic acceleration/deceleration of relativis- 
tic charged particles, thereby affecting the observed inten- 
sity. Change can also occur as a result of the emission region 
passing through a bend in the jet causing, amongst other 
things, a change in viewing angle with respect to the mo- 
tion of the emission region, and hence a change in Doppler 
factor. 

For the examples above, the observer-frame variability 
time depends on the geometry of the emission region (see, 
e.g., Protheroe 2002) and, for viewing the assumed cylin- 
drical emission region along its axis, is at least £' /{2Dc). In 
general, arbitrarily low Doppler factors are possible for ar- 
bitrarily high £' /r. However, in this case this is clearly at 
the expense of having a time scale for variability which may 
be unreasonably large, and of course also requires the jet to 
be extremely well collimated over the length of the emission 
region. 




Figure 6. Observer-frame variability time (thick curves) and 
energy-loss time (thin curves) vs. Doppler factor, D, for: Uc = 
lO^ i, J7 = 1 and (£'/r) = 1 (solid curves); Uc = lO^^'^s, r) = 1 
and {£' /r) = 1 (dotted curves); and Uc = lO"'^'^'^^, r; = 1 and 
{£'/r) = 316 (dashed curves). The variability times for models 
A-O are indicated. 



Alternatively a shock may pass through the emission re- 
gion compressing the magnetic field and accelerating parti- 
cles - a pre-existing mono-energetic population of relativistic 
particles will receive an energy boost by a factor of ~ 7'fhock 
and acquire a power-law tail as a result of diffusive shock ac- 
celeration. In the case of a plane shock moving at jet-frame 
speed /^shockC along the jet (see, e.g. Protheroe 2002), the 
observer-frame variability time in our case {0' = 0) is at 
least 



var, shock 



2c 



^ shoe' 



(46) 



which is longer (shorter) than £' / (2_Dc) if jS'^i-^^ck is less than 
(greater than) 0.5 for a forward-moving shock. 

Except for a point-like emission region following a bent 
or helical trajectory such that the Doppler boosting fac- 
tor, changes rapidly with time, the variability time can not 
be shorter than the observer-frame energy-loss time scale, 
D~^E' /{dE' /dt'). This is easily obtained from equation 1321 
and is plotted (thin curves) against Doppler factor in Fig.|S] 
for the following cases: Vc = 10^'^, 77 = 1 and {£' /r) = 1 
(sohd curves) Vc = 10"■'^^ 7? = 1 and {£' /r) = 1 (dotted 
curves); and Uc = 10"■®^ = 1 and {£'/r) = 316 (dashed 
curves). As noted above, the variability time will in general 
be longer than the energy-loss time scale due to the dimen- 
sions of the emission region. For the present geometry and 
viewing angle, I take the variability time to be 

1/2 

E' 



= D~ 




{dE'/dt') 



(47) 



but one should note that in the case of shock excitation 
that it can be longer than this for /^shock < 0-5, or as short 



as D-^E'/ {dE'/dt') if /3^ho 



1. In Fig. |H] I have also 



plotted tvar vs. D for the same three cases, and the values for 
models A-0 are indicated by the letters. Of the models with 
Doppler factors less than ~ 10'^, models C-E give variability 
times ~ 1-30 years. Models F-J give variability times of ~ 1 
week to ~ month, and are quite compatible with the long- 



10 R. J. Protheroe 



term variability observed for PKS 0405-385. Models K-O, 
with larger £' /r, give variability times of ~ 1 month to ~ 100 
years, with model M having t var ^ 1 ySSJ". 



7 CONCLUSION 

To obtain low Doppler factors, without invoking extreme 
variations from equipartition, requires a larger angular di- 
ameter for the IDV core of PKS 0405-385 than the ~ 
6 micro-arcscc assumed by Kodziora-Chudczer et al. (1997) 
corresponding to a distance of ~ 500 pc to the ionized 
material responsible for the interstellar scintillation. A fur- 
ther reduction in the minimum Doppler factor needed to 
avoid the Compton Catastrophe can be obtained by hav- 
ing optically thin IDV core emission, preferably originating 
from a region elongated along the jet, and observing it at 
a small viewing angle with respect to the jet axis. An an- 
gular diameter a factor of 2 4 larger than assumed by 
Kedziora-Chudczcr ot al. (1997) and Walker (1998) corre- 
sponds roughly to models D, H-I and M-N, all of which 
are allowed by the present data, and have Doppler factors 
~ 200, 80-20 and 30-8, respectively. Interstellar scintilla- 
tion screen distances as small as 20 pc, possibly associated 
with material in the local bubble, have previously been in- 
voked to explain IDV in J1819-I-3845 (Dennett-Thorpe & de 
Bruyn 2000, 2002). By using a distance of 20 pc, even lower 
Doppler factors could easily fit the observations. 

In conclusion, if one observes optically thin emission 
along a long narrow emission region, the average energy den- 
sity in the emission region can be significantly lower than 
(47r/c) / Iiydv. If one takes into account the uncertainty in 
the distance to the ionized clouds responsible for interstellar 
scintillation causing rapid IDV in PKS 0405-385 the bright- 
ness temperature could be £is low as ~ 10^^ K, or lower, 
at 5 GHz. The radio spectrum can be fit by Ii, oc v^^^ ex- 
pected if the emission is optically thin and the electrons are 
mono-energetic, or have a minimum Lorentz factor whose 
critical frequency is well above the frequency range of in- 
terest. Such a spectrum would occur naturally in models 
in which are produced as secondaries of interactions of 
protons, e.g. by Bethe-Hcitlcr pair production, or through 
TT^ decay following collisions of with matter or 

low-energy target photons (pion photoproduction) by pro- 
tons, or by neutrons themselves produced in pion photopro- 
duction interactions (e.g. p7 — »■ nTr"*") perhaps closer to the 
central engine. The combination of all these factors enables 
the Compton catastrophe to be avoided, and also predicts 
that X-ray and gamma-ray emission at observable flux lev- 
els should be expected from compact cores of extreme IDV 
sources such as PKS 0405-385. 
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